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Abstract 

We consider the Schrodinger operator on the real line with a 2 x 2 matrix valued 1- 
periodic potential. The spectrum of this operator is absolutely continuous and consists 
of intervals separated by gaps. We define a Lyapunov function which is analytic on 
a two sheeted Riemann surface. On each sheet, the Lyapunov function has the same 
properties as in the scalar case, but it has branch points, which we call resonances. 
We prove the existence of real as well as non-real resonances for specific potentials. 
We determine the asymptotics of the periodic and anti-periodic spectrum and of the 
resonances at high energy. We show that there exist two type of gaps: 1) stable gaps, 
where the endpoints are periodic and anti-periodic eigenvalues, 2) unstable (resonance) 
gaps, where the endpoints are resonances (i.e., real branch points of the Lyapunov 
function). We also show that periodic and anti-periodic spectrum together determine 
the spectrum of the matrix Hill operator. 



1 Introduction and main results 

We consider the self-adjoint operator Ty = —y" + V(x)y, acting in L 2 (R) © L 2 (M), where V 
is a symmetric 1-periodic 2x2 matrix potential which belongs to the real space T~C p ,p = 1, 2, 
given by 



T V 1 V, 



tc p = \ v = v* = v 1 = ( :; i : / v 3 (x)dx = o 



^3 V 2 J J Q 



equipped with the norm || V||^ = f Q (\Vi(x)\ p + \V2(x)\ p + 2\Vz(x)\ p )dx < oo. Without loss of 
generality we assume 



V{o)= [ V(t)dt = di&g{V w ,V 20 }, V W ^V 20 , V m0 = [ V„ 
Jo Jo 



(x)dx, m = 1,2. 
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Introduce the self-adjoint operator T° = — 4s, with the domain Dom(T°) = W$(M) © 
W|(R). In order to get self-adjointness of T we use symmetric quadratic forms. We 
shortly recall a well known argument (see [RS]). We define the quadratic form (Vip,ipi) = 
— J R Vip fadx, tp,tpi G Dom(T°). Using the estimate (see [Kl]) \(q'f,f)\ < e(f'J) + 
b e (f,f) for any small e > and some b £ > and any / G Q £ L 2 (R/Z) we deduce 

that 

|(V>,Vi)l < (l/2)ftW) + ^ e © W^R). 

Thus we can apply the KLMN Theorem (see [RS]) to define T = — + V. There exists a 
unique self-adjoint operator T with form domain Q(T) = Wf(R) © W^R) and 

(Tip, ipi) = (-rl/'M + (ViPM, all ^1 G Q(T°) = W?(R). 

Any domain of essential self-adjointness for T° is a form core for T. 

It is well known (see [DS] p. 1486-1494, [Ge]) that the spectrum a{T) of T is absolutely 
continuous and consists of non-degenerated intervals S n ,n = 1,2..,. These intervals are sep- 
arated by the gaps G n with lengths \G n \ > 0,n = 1,2.., iV^ ^ 00. Introduce the fundamental 
2 x 2-matrix solutions <p(x, A), $(x, A) of the equation 

-y" + V(x)y = \y, A e C, (1.1) 

with initial conditions </?(0, A) = i?'(0, A) = 0, </?'(0, A) = i?(0, A) = J 2 , where I m ,m > 1 is 
the identity m x m matrix. Here and below we use the notation (') = d/dx. We define the 
4 x 4-monodromy matrix by 

&(x,\) <p(x,\) \ , v 

#>(x,\) <pf(x,\) J ■ [L - Z) 

The matrix valued function M is entire. An eigenvalue of M(A) is called a multiplier, i.e., it 
is a root of the algebraic equation D(t, A) = 0, where 

D(t, A) := det(M(A) - r/ 4 ), r, A G C. (1.3) 

There is an enormous literature on the scalar Hill operator including the inverse spectral 
theory [M], [GT], [KK]. In the recent paper [Kl] one of the authors solved the inverse problem 
(including characterization) for the operator —y" + v'y, where a function v G L 2 (T) and v' is 
the distribution and T = R/Z. However, in spite of the importance of extending these studies 
to vector differential equations, apart from the information given by Lyapunov and Krein (see 
[YS]), until recently nothing essential has been done. The matrix potential poses interesting 
new problems: 1) to construct the Lyapunov function, 2) to define the quasimomentum 
as a conformal mapping, 3) to derive appropriate trace formulae (e.g. analog to the trace 
formulas in [K2]), 4) to obtain a priori estimates of potentials in terms of gap lengths, 5) 
to define and to study the integrated density of states. In fact this is the motivation of our 
paper. 

The basic results in the direct spectral theory for the matrix case were obtained by 
Lyapunov [Ly] and Krein [Kr] (see also Gel'fand and Lidskii [GL]). Below we need the 



M(A) = M(l, A), , M(x,X) 
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following well-known results of Lyapunov [YS], which we formulate only for the case 2x2 
matrices. 

Theorem (Lyapunov ) Let V E H 1 . 

i) If r(A) is a multiplier for some A G C (or X E WL), then t _1 (A) (or r(A) ) is a multiplier 
too. 

ii) M(A) has exactly four multipliers Ti(A), rf 1 (A), r 2 (A), t 2 _1 (A) /or a// A E C. Moreover, 
A G o-(T) z/f |r m (A)| = 1 /or some m = 1, 2. 

jhJ // r(A) zs a simple multiplier and |r(A)| = 1, i/ien r'(A) 7^ 0. 

We mention the papers relevant in our context, for various problems periodic systems we 
refer the reader to [CGHL]. In the paper [Cal,2] Carlson obtained the trace formulas. In 
[Ca3] he proved the compactness of Floquet isospectral sets for the matrix Hill operator. 

By the Lyapunov Theorem, each M(A), A G C has exactly four multipliers r m (A), r m 1 (A), 
m — 1,2, which are the roots of the characteristic polynomial det(M(A) — r/ 4 ) = 0. If 
c = y20 2 Vl " > 0, then the multipliers have the following asymptotics 

r m (A) = e^-^ + ° (1/A)) , IAI-00, \V\-*n\>^m = l,2, (1.4) 
(see Lemma 3.2). Next we define the functions 

Note that <p(x, A), A), pi(A), P2(A), p(A) are real for all x, A G 1R and entire. If c = 
v ' 20 ~ Vw > 0, then Lemma 3.1 yields p(A) = pg(A)(l + o(l)), as |A| -> 00 in X>i, where 



: = |A G C : |A| > r 2 , (v^ - im\ > j,n E n|, r > 0. (1.6) 



Thus we may define the analytic function y/p(X), A G XV, for large r > 0, by the requirement 
that a/ p{t) = po(A)(l + o(l)) as |A| — > 00 in D r . Then there exists a unique analytic 
continuation of ^/p from V r into the two sheeted Riemann surface A (in general, of infinite 
genus) defined by yfp. We now introduce our Lyapunov function A (A) by 

A(A)=/i 1 (A) + v / pTA), AG A. (1.7) 

Let Ax(A) = pi(A) + \j 1 p{X) on the first sheet Ai and let A 2 (A) = pi(A) — y/ p(X) on the 
second sheet A 2 of A. Now we formulate our main result concerning the function A. 
Theorem 1.1. Let V E Ti. 1 with V10 < ^20 ■ Then the function A = p\ + ^fp is analytic on 
the two sheeted Riemann surface A and has the following properties: 

i) 

A m (A) = rm(A)+ 2 Tml(A) , A G A m , m = l,2, (1.8) 

sin \/~X /p\^ m ^\\ 
A m (A) = cosv / A + K„o^^ + o(^^J, m = l,2, A G V 1 . (1.9) 
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ii) X G C belongs to cr(T) iff A m (X) E [—1, 1] /or some m = 1,2. 
zzzj If X E a{T), then p(X) ^ 0. 

("ZTie monotonicity property.) Let A m be real analytic on some interval I = (cti, ct 2 ) C M 
and —1 < A m (A) < 1, for any X E I for some m E {1, 2}. T/ien A' m (A) 7^ /or eac/i X E I . 
Remark, i) For the scalar Hill operator the monodromy matrix has exactly two eigenvalues 
r and r -1 . The Lyapunov function |(r + r _1 ) is an entire function of the spectral parameter 
and it defines the band-gap structure of the spectrum. By Theorem 1.1, our Lyapunov 
function for the matrix Hill operator also defines the band-gap structure of the spectrum, 
but it is the sheeted analytic function. 

ii) Consider the case of a diagonal potential, i.e. V3 = 0. Then the Riemann surface 
degenerates into two components, and we get 

= ^(A(D + A( 2) ), P = ^(A(i) - A (2) ) 2 , 4~p = -(A (1) - A (2) ), (1.10) 

where A( m ) is the Lyapunov function for the scalar Hill operator — y" + V m y, m = 1, 2. Thus 
Ai = Pi + y/p = A(i) and A 2 = \i X - y/p = A (2) . 

iii) A lot of papers are devoted to the resonances for the Schrodinger operator with 
compactly supported potentials on the real line, see [BKW],[K3],[S],[Z]. Assume that we have 
the couple constant before the potential. If this constant changes then roughly speaking 
some resonances create the eigenvalues. In our case, if the couple constant (before the 
periodic matrix potential) changes, then roughly speaking some resonances create the gaps, 
see Proposition 1.3. 

Let D±(X) = jD(±l, A). The zeros of D + (X) and D_(X) are the eigenvalues of periodic 
and anti-periodic problem associated to the equation — y" + Vy = Xy. Denote by A 2n fc,n = 
0,1,.., and k E {1,2,3,4} the sequence of zeros of D + (counted with multiplicity) such 
that A 0) i ^ A , 2 ^ A 2j i ^ A 2j2 ^ A 2j3 ^ A 2>4 ^ A 4) i ^ A 4j2 ^ A 4j3 ^ A 4j4 .... Denote by 
A2n-i,fe> {. n i k) E N x {1, 2, 3, 4} the sequence of zeros of _D_ (counted with multiplicity) such 
that Ai,! ^ Ai )2 ^ Ai ; 3 ^ Ai j4 ^ A3,i ^ A3 j2 ^ A3 t 3 ^ A3 j4 . . .. Note that A^fc, n — 0, 1, .., and 
k E {1, 2, 3, 4} are the eigenvalues of problem with period 2 for the equation — y" + Vy = Xy. 

Denote by {r n }f the sequence of zeros of p in C (counted with multiplicity) such that 
^ |ri| ^ |r 2 | ^ |r 3 | ^ . . .. We call these zeros of p the resonances of T. We formulate the 
theorem about the recovering the spectrum of T and the asymptotics of the periodic and anti- 
periodic eigenvalues and resonances at high energy. Furthermore, we write a n = b n + £ 2 (n) 
iff the sequence {a n — b n } n ^i E t 2 . Recall that Vo, V m0 , m = 1, 2 are defined by (jl.5j) . 
Theorem 1.2. Let V E H 2 with V20 — Viq > 0. Then the following asymptotics are fulfilled: 

X n ,m+k = (t™) 2 + V m o + f (n) , m,k = l,2, U -> OO, (1.11) 

r2n-m = (nn) 2 + £ 2 (n), m — 0, 1, n — > oo. (1-12) 
ii) Let V E Ii 2 , V 20 — Vio > 0. Then the following statements hold 

a) The periodic spectrum and the anti-periodic spectrum recover the resonances and the 
spectrum of the operator T . 

b) The periodic (anti-periodic) spectrum is recovered by the anti-periodic (periodic) spectrum 
and resonances. 
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Less precise asymptotics for the case V G C 2 were obtained in [Cal]. 
Example. Consider the operator T ljV = — + g 7i ^, where g 7il/ = a J + 7t>„(x)J 1 is a 

potential for some constants a, 7 G K, 

J= C 1 ,),J 1= (° M , and , w satisfies 



x -1 y ' 1 V 1 , 
Condition A. ^ — n a G (0, 1) /or some integer n a ^ 0. Each function v v G L 1 (T), v = 
1, zs suc/i i/iai = iv(l — x), x G (0,1), and J" v v {x)dx = 1 and for any 

f G C(0,k),k G N ; i/ie following convergence holds: 



r-k rk ^ j 

/ v v {x)f{x)dx — > / 5 per (x)f(x)dx as 1/ — >■ 0, (5 per = 2> <5(x — n — -). (1-13) 

If 7 = 0, then we have the operator T° = — y" + g° with a constant potential q° = aJ. 
In this case there are no gaps in the spectrum. The fundamental solutions of Eq. fjl.lj) with 
q° = a J have the form $° = diag(c + , c_), (p° = diag(s+, s_), where 

c±(x,X) = cosi]±x, s±(x,X) = - — r]± = \/\ =F a, diag(6, c) = ( J " ] 

*7± V u c / 

and the branch of a/A is given by vT = 1- Below we sometimes will write p(A, V), M(A, V), .., 
instead of p(A),M(A), .., when several potentials are being dealt with. Then the functions 
P°(') = P('i A*m(') = / i m(',Q , °), ■• corresponding to g° have the forms 

o m c + (m,A) + c_(m,A) (c+(l, A) - c_(l, A)) 2 

L>i(A) = (lTc + (l,A))(lTc-(l,A)), A° 1} (A) = c+(l, A), A° 2) (A) = c_(l, A). (1.15) 

Here A9 x is the Lyapunov function for the Eq. — y" — (— l) m ay = Xy, m = 1, 2. The entire 
function c + (l, A) — c_(l, A) has only the simple roots z® given by 

z° n = (Tm) 2 + -^ = Ta + (nn±^-) 2 , n>\. (1.16) 

Note that z® > z\ > ... > z® a and z® a+1 < z® a+2 < .... Thus all roots of p° = p(-,q°) 
have multiplicity 2 and are given by (|1.16j) . If |§ £ N, then the zeros of A9,(A) 2 — 1 have 
multiplicity 2 and are given by 

A^ >n = a - (-ir(Trn) 2 , A° n ^ A° s , m=l,2, n, O 0. (1.17) 

Note that A^ n , to = 1, 2, n ^ are the roots of D±(X). If 7, ^ 7^ are small enough, then 
there exist gaps in the spectrum of T 7j ^. Define the disk D r = {A : |A| < iT 2 r 2 },r > 0. We 
show that there exist the non-degenerated resonance gaps for some V. In this example some 
resonances r n G R and some r n ^ M. 

Proposition 1.3. Let a potential q liV = aJ + jv v (x) J\, a > 0, 7 G M, satisfy Condition 
A, and let ^§ ^ N. T/ien /or any ^ 1 + a t/iere ea;zs£ small e, 61,62 > snc/i t/iat 
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for any (<y,v) G (0,£i) x (0,e 2 ) all zeros z±(q 7:U ) of p(X,q 1)U ) and the zeros A* )Tl (g 7)1/ ) of 
A^(A, q lih ) — 1 in the disk D/v+i are simple and have the properties 

z^ w )6C ±1 Un<n„, z~(g 7i „) < «+(g 7)I/ ), n a < n < iV, |^(g 7 , v ) - z°\ < e 2 , 

(1.18) 

A m,n(?7,^) < A m,n(?7,^ l^nfo.") ~ A m,nl < 6 2, Tl = 0, .., N, m = 1, 2. (1.19) 

T/iere are no o£/ier roots of p(\, g 7) „) and A^(A, q ltV ) — 1 in the disk B> N+ i. 

Remark, i) If < a < 2ir 2 , then p(X,q J U ) has only real roots zf(q ltV ) in the disk B> N+ i for 

small 7, za If a > 2n 2 , then p(A, g 7ii/ ) has at least two non-real roots z x {q ltV ) for small 7, za 
ii) We show that the operator T 7jJ/ has the gaps associated with the periodic or anti- 
periodic spectrum. Moreover, we show the existence of new gaps (so-called resonance gaps). 
The endpoints of the resonance gap are the branch points of the Lyapunov function, and, 
in general, they are not the periodic (or anti-periodic) eigenvalues. These endpoints are 
not stable. If they are real (n = n a + 1, ..,N), then we have a gap. If they are complex 
(1 ^ n ^ n a ), then we have not a gap, we have only the branch points of the Lyapunov 
function in the complex plane. 



2 Fundamental solutions 

In this section we study 1?, tp. We begin with some notational convention. A vector h = 
{hn}\ £ has the Euclidean norm \h\ 2 = Yli \^ri\ 2 , while a. N x N matrix A has the 
operator norm given by \A\ = sup| h | =1 \ Ah\. The function tp satisfies the following integral 
equations 

/ f x sin vX(x — t) Tr , s , , s , , , s sinV^r: ,„ „. 

p(x,X) = <p Q (x,\) + J -£= J -V(t)tp(t,X)dt, tp (x, X) = — -j= — J 2 , (2.1) 

where (x, X) G R x C. The standard iterations in (|2.1|) yields 

<p(x, X) = J2 <Pnfa X) , p n+1 (x } X) = [ X %~^ V(t)<p n (t, X)dt. (2.2) 

The similar expansion -& = ^2 n> Q& n with $ (x,X) = (cosvAa;).^ holds. We introduce the 
functions 



i£(A) = / dt cosVX(m-2t + 2s)F(t,s)ds, F(t, s) = Tr V(t)V(s), (2.3) 
Jo Jo 

m = 1, 2. In Lemma 2.1 we shall show the simple identity 

J 2 °(A) = 41° (A) cos y/X. (2.4) 

We define |A|i = max{l, |A|} and 



V (0 )= [ V(x)dx, V (1) = TrV (0) , V (2) = Tr V 2 Q) , A = e l Im ^ + ", x=^=. 
Jo 



\V\\i 



We prove 

Lemma 2.1. For each (x, V) G K + x Ti 1 the functions (p(x, -),i)(x, ■) are entire and for any 
N ^ — 1 the following estimates are fulfilled: 



max 



N N N 

{ \#(x, A) - M*, A) | , \V\{ip(x, ^(x, A)) | , | -j= (V(x, A) - J] A) 



N 



N+l 



(N + l)\ 



A x . 



(2.6) 



Moreover, each fi m {X),m — 1,2 is real for A GR and entire and the following estimates are 
fulfilled: 



\li m (X)\^A m , 

sin mz (rax' 
|/^m(A) — cosmz ; raVm I ^ — ; — A 



A^m(A) — cosraz| ^ raxv4 r 

I / / / x ; " 

(1)1 



4z 



|/^m(A) — cosraz 



sinmz Tr 1 , Tn .,, m'cosmz Tr ,. (rax) , ,, 

-^(1) - ^( J ™( A ) n ^(2))l < y -^A m , 



where 7°, 7° satisfy (12. 4j) and z = Va. 

Proof. We prove the estimates of the proof for ip', is similar. (|2.2jl gives 



V9 n (x, A) 



<po{xi,x)[ n 



sin(\/A(x fc+ i - x k )). 



0<a;i<a;2<---<a:n+l=a; 



Va 



(2.7) 
(2.8) 

(2.9) 



)dxidx2---dx n , (2.10) 



where for matrices 01,02, ...,an we denote \\ a^ = a n a n -\...a\. Substituting the estimate 

|VAyo(«, A) I < el 1 " 1 ^^ into (l2~TUJ) we obtain |VA^(x, A)| sC el Im2 l x , which shows that 
for each x ^ the series (J2.2J) converges uniformly on bounded subsets of C. Each term of 
this series is an entire function. Hence the sum is an entire function. Summing the majorants 
we obtain estimates (|2.6j) . 

We have 4// m = Tr M m (A) = Tr M(ra, A) = Tr J2 n ^o M nijn, A), where ra = 1, 2 and 

Tr M (m, A) =4cosmz, Tr M n (m, A) = Tr$ n (ra, A) + Tr<// n (ra, A), n > 1. (2.11) 
The estimates \tp' n {m,\)\ < ln^ e l Im2 l m and |tf n (ra,A)| ^ (n^ e l Imz l m yield 



Tr M n (m, A) | ^ 4 



,mx 



n! 



3 m| Im %/A| 



, n ^ 0. 



(2.12) 



Using (|2.11|) we obtain 



Tr Mi(m, A) 



sinz(ra— t) coszt+cos,2(ra— £) sinzt) Tr V(t)dt = ml^i), (2.13) 



7 



and 



1 

since 



o Jo 



j rm rt 

Tr M 2 (m, A) = — / sin z(t — s)z(m — t + s)F(t, s)dtds 
z Jo Jo 

/ (cos z(m — 2t + 2s) — cos zm)F(t, s)dtds = — ^ (i^(A) — cos mzV^ 

Jo 2z 2 



J F(t,s)dtds =2 Tt \J Vi $ dt ) = T^ (2) ' 



'o jo 

We have that fix, /i 2 are entire. Moreover, the trace of the monodromy matrix is a sum of its 
eigenvalues. By Lyapunov Theorem (see Sect.l), the set of these eigenvalues is symmetric 
with respect to the real A £ 1. Hence fix, fi 2 are real on R. 

We show (HH). Let J m = J mO)- We have 2 cos zlf =Y + Yl, where 

y = f f cos2z{t - s)F(t,s)dtds, Y x = [ [ cos 2z(l - t + s)F(t, s)dtds. 
Jo Jo Jo Jo 

We get 1$ = Yx + Y 2 + Y 3 , where 

Y 2 = f f cos 2z(l -t + s)F(t,s)dtds, Y 3 = I f cos 2z(l - t + s)F(t, s)dtds 
Ji Jo Ji Ji 

and using the new variable r = t— 1 we get Y 2 = JJ" cos 2z(r — s)F(r, s)drds = 2Y . We use 
the new variables r = t—1, a = s — 1 and obtain Y% = f J Q T cos2z(l—r+cr)F(r, a)drda = Yx, 
which yields 1$ = 2Y X + 2Y . Thus we have (J23J). ■ 

We need the basic properties of the monodromy matrix. 
Lemma 2.2 Let V G H 1 . Then the function D(t, A) = det(M(A) - r/ 4 ), A, r 6 C 2 is entire 
on C 2 and the following identities are fulfilled: 

D' T (r,X) = -D(r,X)Tr^M(X)-rI^ \ (2.14) 

D(t, ■) = r 4 - 4/iir 3 + 2(4/i 2 - /i 2 )r 2 - Afi X T + 1 
= (r 2 - 2(fix - Jp)T + l) (r 2 - 2(/i! + y/p)r + l) . (2.15) 

Proof. Let D(r) = D(r, A). The standard identity D'(t) = D(t) Tr((M -tT 1 ^^) yields 
(l2~T4l . We prove (l2~T5T) . Since det M = D(0) = 1, we have: D(r) = 1 + ar + br 2 + cr 3 + 
<ir 4 , a, b,c,d G C. Then -D(r) = (r — 7~i)(t — t 1 _1 )(t — r 2 )(r — r 2 _1 ), where r 1; rf , r 2 , r 2 _1 are 
the multipliers. Therefore d = 1, a = c. Then we have 

£>(t) = 1 + D'(0)t + ^L>"(0)r 2 + D'{0)t 3 + r 4 . (2.16) 

Identity (pHIjl yields D'(0) = - Tr M _1 = - Tr M = -Afi x . Differentiating (j2~T4jl we obtain 
D"(0) = -D'(t) Tr(M - r)" 1 - D(r) Tr(M - r)- 2 | r=0 = 4(4/i 2 - /i 2 ). Substituting these 
identities into (|2.1fi|) we obtain the first identity in (|2.15jl . The second identity is proved by 
direct calculation. ■ 
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3 The Lyapunov function 

We need some results about the functions p, p = c o^0^> where Co = |(^20 — ^io)- 
Lemma 3.1. i) For each V G Ti 1 the function p = + 1) — p 2 is entire and real on the 
real line. Moreover, the following estimate is fulfilled: 

|p(A)-pg(A)| ^2x 3 e 2 l Im ^ +2 -, AgC, x=||y|| 1 /vlAl7. (3.1) 

ii) Let V G Ti. 1 with Cq = |(V2o — Vio) > 0. Then for each integer N > 2 ^J^ 1 the function 



p(A) has exactly 2N roots, counted with multiplicity, in the disk {A : |A| < 7r 2 (iV+|) 2 } and for 
each n > N , exactly two roots, counted with multiplicity, in the domain {A : |vA — 7m| < 1}. 
There are no other roots. 

Hi) Let V G W x ,co > 0. Then the function Jp~ is an analytic function in the domain 

2 9 |[V|| 3 

T> r , r = — given by (jl.fij) and the following estimate is fulfilled 



|VpTA)- Po (A)|<^^#, C ^4 4 ^l<Vp, XeVr _ (3 _ 2) 



5 x/|A[' q 2 



o * 

Proof i) By Lemma 2.2, p is entire and real on the real line. Let p m = p m o+/ i mi+Pm2+/^m3, 
m = 1, 2, where 

Sill TflZ 1 77?^ 

p m0 = cosm2;, p m i = 4 ^ mV(i), p m2 = ^(/^(A) - — cosmzV (2 )), (3.3) 

and z = y/X, where V( m ) is defined by (12.5)) . We obtain p 2 = £>i + £> 2 , where 

Si = (pio + AHi) 2 + 2pi pi2, ^2 = 2^ioPi 3 + (/ii - + pn)(pi - pio - Pn)- 
Then (13.3)1 yields 

sin 2z 4 cos zi^ (A) + sin 2 zV 2 } - 2 cos 2 zV {2 ) 



ft = cos' , + + - 16 /' . 

Thus we get 

^2 + 1 2 ^ . ^ ~ 1 + ^20 + P21 + P22 „ n P23 p 

P = — ^ /U x = Cri + (jt 2 , Cri = £>i, G 2 = — i>: 



2- 



J 2 °(A) - 2 cos 2zV(2) - 4 cos z/°(A) + 2 cos 2 zV {2) - sin 2 zV ( 2 } s i n 2 z 2 

Gl = To? = i6^ (2y(2) " 

which yields G x = p§. Using ()2.7)) - (|2.9)) we obtain 

|52|^2^ 2 + (xA + xA)^ 2 = ^ 2 , 
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\G 2 \ ^ ^ + \B 2 \ ^ ^-A 2 + —A 2 = 2x 3 A 2 , (3.4) 
2 3! 3 v ' 

which yields ()3.1|) . 

ii) Let N' > N be another integer and r = tc(N + |). Note that x ^ sTf^jjVjp ^ 2^ ^ or 
A G P r . Using the estimate e' Imv/ ^l < 4| sin \/A|, AeD r and (|3.1|) we obtain (on all contours) 

\p(X)-pl(X)\^2^e 2 ^^^^\pl(X% C = 4^<^±, \eV r . (3.5) 







A 



Hence, by the Rouche theorem, p(X) has as many roots, counted with multiplicity, as 

in the bounded domain V r fl {A : |A| ^ tt(N' + §)}. Since sm2 , has exactly one double root 
at (7m) 2 , n ^ 1, and since N' > N can be chosen arbitrarily large, ii) follows, 
iii) Let p = p 2 , + Pi, Po = c o S f^- Estimates ()3.5j) imply 

V^A) = Po(A)v / T+&(A), 6 = ^, |6(A)|<^<J, A G V r . (3.6) 



Using the estimate | ^/l + 6(A) — 1| ^ §|fr(A)| f° r A G Z> r , we obtain 

IVKA) -Po(A)| = |p (A)(v / rT6(A) - 1)1 < ■ 

We need asymptotics of the eigenvalues of the monodromy matrix. 
Lemma 3.2. Let V G Tl 1 , Vq > 0. T/jen the monodromy matrix M(A) /ias two multipliers 
T m (X),m= 1,2 wif/i asymptotics (jl.4|) . Moreover, asymptotics ()1.9|) is fulfilled. 
Proof. Using estimates (|2.8j) . (|3.2j) we obtain asymptotics (jl.9j) . Asymptotics (jl.9j) yield 

A^(A) - 1 = (cos V\ + V m o S -^) 2 - 1 + 0(E(X)) 

9 j— sin V'A cos v'A 9 /—( Tr cos \[X „ . „ , , . . 

= - sin 2 >/X + Kno rr + °(^( A )) = " sin ^ 1 " y mo r- ^ + 0(E(X)) 

V A v v A sin v A 

as A G Pi, |A| — > 00, where £>(A) = ^''T^ 1 , which implies 



V^A) - 1 = i sin v^A - iVmo^P + 0(E(X)). (3.7) 

2v A 

By ()2.15p . the matrix M(A),A G V r , has the eigenvalues r m (A) satisfying the identities 
r TO (A) + r m (A) _1 = 2A m (A). Then r m (A) has the form r m (X) = A m (A) + a/A^A) - 1 and 
the asymptotics give 

r m (A) = + K»of^5 - + OW)) = e^(l - ^)+0(£(A)) 
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which yields (jl.4j) . ■ 

Now we prove our first result about the Lyapunov function A = /ii + ^fp. 
Proof of Theorem 1.1. From Lemma 3.1 we obtain the analytic function A on the 
Riemann surface of the function ^fp. 

i) Identity (J2.15J) shows that A m = for some multiplier r. Lemma 3.2 gives the 
asymptotics of A m and r m , m = 1,2. 

ii) The result follows from the statement i) and the Lyapunov Theorem (see Sect.l). 

hi) If A G cr(T), then pi(A) is real. By ii), A(A) is also real. Hence by (jl.7j) . a/ p(A) is real 
and p(A) ^ 0. 

iv) Assume that A' m (Ao) = for some Ao G / C o~(T) and m G {1, 2}. Then we have 

a(p)/\ \ 

A m (A) = A m (A ) + (A - A ) p - \ 0) + 0(\X - X \ p+1 ), as A - A -> 0, (3.8) 

pi 

where A#(A ) ^ for some p ^ 1. By the Implicit Function Theorem, there exists some 
curve r C {A : |A — Ao| < e} fl C+, r ^ 0, for some e > such that A m (A) G (—1, 1) for any 
A G T. Thus we have a contradiction with the Lyapunov Theorem in Sect. 1. ■ 

Recall that D±(\) = ^det(M(A) =F ^4), the set {A : -D+(A) = 0} is a periodic spectrum 
and the set {A : -D-(A) = 0} is an anti-periodic spectrum. Now we prove a lemma about the 
number of periodic and anti-periodic eigenvalues in a large disc. 

Lemma 3.3. For each V G Ti 1 the functions Ai + A 2 , AiA 2 , D± are entire and satisfy the 
following identities: 

A? + A^ = 1+^2, AiAj, = 2/i? - (3.9) 
D ± = ^ lT lf- P = ( 2 ^1) 2 -/^ D+ _ D= _ A ^ (310) 
Let in addition Tr V{t)dt = 0. Then the following estimates and properties are fulfilled: 

maxj (A) -4 sin 4 ^5 , £L(A)-4 cos 4 ^ } ^ x 2 (2 + x) V 1 a £ c (3 n ) 

i) For each integer N > 8||V||i the function D + has exactly 4iV + 2 roots, counted with 
multiplicity, in the disc {|A| < 47r 2 (iV-|- ^) 2 } and for each n > N , exactly four roots, counted 
with multiplicity, in the domain {\\/~\ — 2im\ < f }■ There are no other roots. 

ii) For each integer N > 8j| V^|| 1 the function D_ has exactly 4N roots, counted with mul- 
tiplicity, in the disc {|A| < Att 2 N 2 } and for each n > N, exactly four roots, counted with 
multiplicity, in the domain {\V\ — ir(2n + 1)| < f }. There are no other roots. 

Proof. By Lemma 2.1, 2.2 the functions Ai + A 2 , AiA 2 ,.D-t are entire and identities ()3.9|) . 
(13.10)1 are fulfilled. Using Lemmas 2.1 and 3.1 we obtain 

2 2 
pi(\) = cosz + pi 2 (A), |/Z 12 (A)| < y A2 ' l^ A )l ^ ( 1 + 2x )y A2 - 
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Substituting these relations into D± = (/^ =)= l) 2 — p we get (j3.11l) . 

i) Let N' > N be another integer. Let A belong to the contours C (2N + 1),C (2N' + 
l))C2n(|)j \ n \ > N, where C n (r) = {A/A : |vA — im\ = Tcr},r > 0. Note that x ^ on all 
contours. Then (j3.11|) and the estimate e2 Imv/ ^ < 4| sin ^| on all contours yield 

n r\\ a ■ 4^ e 2 "(2 + x) 2 2{lmVX{ 1 . 4 VA 

D,(A)-4sm 4 ^ ^— r-^-e 1 v 1 < - 4 sin 4 

+ v ; 2 4 4 tt 2 4 2 

Hence, by Rouche's theorem, D + (\) has as many roots, counted with multiplicities, as sin 4 ^ 

in each of the bounded domains and the remaining unbounded domain. Since sin 4 ^ has 
exactly one root of the multiplicity four at (27m) 2 , and since N' > N can be chosen arbitrarily 
large, i) follows. The proof for Z}_ is similar. ■ 

We are ready to prove Theorem 1.2. 
Proof of Theorem 1.2. i) It is enough to consider the case Vio = cq = — V^o- We prove 
asymptotics (jl.llj) for \2n,m, 1 ^ m ^ 4. The proof for A2 n +i,m is similar. Firstly, we prove 
the rough asymptotics of A2ri., m ,r n . Lemma 3.3.i yields \J\2n,m = 27m + e n , \e n \ < 1, as 
n — > 00. By Lemma 3.3, -D+(A) = 4sin 4 ^ + 0{n~ 2 ) as \yf\~ — im\ ^ 1, n — > 00. Then 
identity -D+(A 2n ,m) = implies 

y/X 2n ,m = 2™ + e n , e n = 0(n- 1/2 ) } 1 ^ m < 4. (3.12) 

Lemma 3.1.ii implies A /r 2n _ m = rm + 5 n , \S n \ < 1 for n —>■ 00, m = 0, 1. Moreover, Lemma 
3.1.i gives p(A) = Cq sm l^ + 0(r?r 3 ), |\/A — 7m\ < 1 as n — >■ 00. Since p(r n ) = 0, we have 

^/r 2n -m = 7m + S n , 5 n = 0{rT 1/2 ), m = 0, 1. (3.13) 

Secondly, in order to improve these asymptotics of \2n,m,r n we need asymptotics of the 
multipliers in a neighborhood of the points 7m. We introduce the matrix M = U~ X MU with 

the same eigenvalues, where U = y J j . We shall show the asymptotics 

M(A) = M (A) + £i(n), M (A) = f f^) £(A) ) ' ^ = ™ + °(^ V2 )' ( 3 - 14 ) 
where 



C(A) = diag(cos?7 + , cosr/_), •S'(A) = diag(sin^ + , sin^_), r/ ± = a/A =F c . (3.15) 
Estimate (|2.6|) gives 

1 f 1 

#(1, A) = cos H 7= / sin V^(l - t) cos VXtV(t)dt + 0(r?r 2 ) 

V A Jo 

\/A + sin v r \{l-2t))V{t)dt + 0{n- 2 ), \V\-7m\^l, (3.16) 
2v A Jo ^ ' 
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cos 



as n — > +00. Let \/\ = nn + u n , u n = 0(n -1 / 2 ). The Taylor formula gives sin2t(7m + u n ) = 
sin27mi + 2tu n cos2imt + 0{n~ 1 ) and the similar formula for cos2t(7m + u n ). Substituting 
these asymptotics into (|3.16j) we obtain 



sin y/X 



tf(l, A) = cos + ^c J + Z\{n) = C(A) + ^(n), = ttw + O^ 1 / 2 ). 

2v A 

Similar arguments for ip(l, A), A), </?'(!, A) yield 

^(l,A) = -^S(A)+£ 2 (n), ^ / (l,A) = -v / A5(A)+£ 2 (n), <//(l, A) = C(A) + £j(n), 

as \/A = 7m + 0(n -1 / 2 ). Substituting the obtained asymptotics into the definition (jl.2|) of 
M and using the identity M = U^MU we get (ETTIJ) . 

We will use standard arguments from the perturbation theory (see [Ka,p.291]). Let A, B 
be bounded operators, A be a normal operator and a(A),a(B) be spectra of A, B. Then 
dist{cr(A), cr(.B)} ^ — Note that M (A) is a normal operator having eigenvalues e tv± . 
Hence M has eigenvalues t± satisfying the estimates \t± — e tr,± \ < \M — M \. Then f!3.14|) 
implies 

r± (A) = e iv± +£1(71), as V\ = nn + 0(r^ 1/2 ). (3.17) 

Now we improve asymptotics (|3.12|) for A 2n/m . Note that A = A 2ni m iff r +(A) = 1 or 
r_(A) = 1. Then (l3~T7l) and (l3~T2D yield 

e ^ ± (A 2 „, m) = x + £2(n); A = x 2nm . (3.18) 
Substituting asymptotics (|3.12j) into (|3.18j) we have 



a/ (27m + e n ) 2 ± c = 27m + ^ 2 (n), n 



-00. 



and therefore e n = =F^; + ^i(w). Substituting this asymptotics into ()3.12|) we obtain p. Ill) 

for \2n,m- 

We improve asymptotics (J3.13j) of r n . Note that A = r n , iff the following condition 
is fulfilled: t+(A) = rZ (A) or r+(A) = t_(A). Asymptotics (|3.17|) imply that the second 
equation has no solutions for large n. We rewrite the first equation in the form: 

j V+ = e - in - + £ 2( n ^ ag y/x = m i + 0{n~ 1/2 ). (3.19) 
Substituting asymptotics ()3.13|) into ()3.19|) we obtain 

y/ [nn + S n ) 2 - c + a/ (7m + 5 n ) 2 + c = 27m + ^(n), 

which gives 5 n = £\{n). Substituting this asymptotics into ()3.13|) we obtain ()1.12j) . 
ii) a) Assume that we have the periodic spectrum Ao, m , m = 1,2, A2 n ,m, m = 1, 2, 3, 4, n ^ 1. 
Using the asymptotics (jl.lljl and repeating the standard arguments (see [PT,pp. 39-40]) we 
obtain the Hadamard factorization for the function D + : 

n 1 n _ (A — Ao,i)(A — A 0) 2j T r A2n,m — A 
D+{X) - 4 11 (27m) 2 • 

m=l,2,3,4,n^l v ' 
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In a similar way, we determine D_ by the anti-periodic spectrum. Using (j3.10j) we obtain p. 
Thus, we recover the resonances. 

b) Suppose, that we have the periodic spectrum and the set of the resonances. Then 
we determine the functions D + by the periodic spectrum and p by the resonances. Using 
(I3.10|) we get pi, P2 and then Z}_. Thus, we recover the anti-periodic spectrum. The proof 
of another case is similar. ■ 



4 Example 



1. Periodic 5-potentials. Consider the operator T 7 = — + g 7 ,7 G C, where q 1 = 
aJ + •ySperJi and the potential S per = X^oo ~ n ~ \)- Let p\ = pi(-, g 7 ), p 7 = p(-, g 7 ), ... 
Lemma 4.1. For the operatorT 7 = — ^ + aJ + ^5 per Ji the following identities are fulfilled: 

p[ = pl pl = pl + 2h, p^ = p° + h, Dl = D° ± -h, h = ^s + s_, 7 eC. (4.1) 

Proof. The solution y(x), x G R of the system — y" + q 1 y = Xy is continuous and y'{x n + 
0) — y'{x n — 0) = •yJiy(x n ) at the points x — x n — n + |. Then the fundamental solutions 
have the form: ^ 7 (x) = i? (x), <p 7 (x) = <p (x), ^ x < |, and 

1 ; V 7c+(i)s-(x-|) c_(z) 7 2 ^ 2' 

V 7s+(|)s-(^- i) s -0) y ' 2 2' 

and 

c+(x) + 7 2 c + (i)s_(l) S+ (x- |) 



i? 7 (x) 



c_(x)+ 7 2 C-(|K(l) S _(x-|; 



s+(z) +7 2 s + (|)s_(l)s + (x- |) * 



for | ^ x < |, where * is some term. These relations yield 

M C+<1) Mi))' ^»-( C+<1) 4) 

M c+(2)+2 \_ (2 ; +2 ,). ^-( c+(2)+2ft M 2 ; +2 , 

The last identities and ()1.14j) imply 

_ Tr(l 7 (l) + (y? 7 )'(l)) _ _ Tr(^ 7 (2) + (y 7 )'(2)) _ 
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which give p 7 = |(1 + 14) - {p[) 2 = p° + h and Dl = (fij T I) 2 - ^ 2 = D° ± - h. * 

We describe the spectrum of the operator T 7 . 
Lemma 4.2. Let the operator T 7 = — + g 7 ; where q 1 = a J + ^8 per Ji, a > and 7 £ R. 

ij Lei < 2^ — n a < 1; /or some integer n a ^ 0. Tnen /or eacn n £ N tnere exist analytic 
functions z^(j), \ j\ < 7 n /or some 7„ > sncn £/ia£ 2^(7) is ine .zero o/p 7 (A) and 

4(7) = 4 ± (-l) n 7^"(c„ + 0( 7 )), c n > 0, k n = I J I ^ ^ na ' , 7 - 0. (4.2) 

Moreover, each spectral interval (z~ (7), z£ (7)) C 1, n > n a is o gap in £ne spectrum ofT 1 . 
ii) If % ^ H, £nen /or eaca n ^ 0,m = 1,2 iaere exisi rea/ analytic functions A^j n (7),7 £ 
(— 7„,7n) /or some 7„ > snca taat A,^ m (7) is the zero of the function A, 2 n (A) — 1 and 

A «, m (7) < A+ m (7), A± m (0) = A° m . (4.3) 

Moreover, each spectral interval (A^ m ( 7 ), A+ m ( 7 )) 7^ 0,n ^ 1 has multiplicity 2. 
Proof, i) Recall that n± = VA =F a. The zero of p 7 (A) = satisfies the equation 

= P 7 (A) = / 2 (A) + 2_s + (A)s_(A), / = 2 sin *±^IL sin 2t±*i (4.4) 
The zeros z® of p° = p(-, g°) have the form (Jl.lfcij) and satisfy the following identities 



V / ^° + a+ v / - 2 n _a = ~; A/z° + a- A/,2;°-a = 27rn, if 1 < n < n a , (4.5) 



V / ^ + a+ v / ^°-a = 27rn, v / ,2°+a- v / ,2°-a = — , if n > n a ^ 0. (4.6) 



7rn 



Using (|4.5jl . (jUBj) we have the identity 



. 1 — cos — > if 1 ^ n ^ n a , 

sin77+(A) sinr/_(A) U=^o = { a nr \ ^ n •* ^ . 4.7 

m-v y / v /ia « n ^ cos — — 1 < if n > n n v ' 



irn 



Recall that the function / has only simple zeros A = n ^ 1. Consider the case n > n a , 
the proof for 1 ^ n ^ n a is similar. We rewrite p 7 (A) = in the form 



(/(A) - 7 ^(A))(/(A) + 7 F(A)) = 0, F(A) = v / - S+ (A)s_(A)/2, F(^) > 0. 

Applying the Implicit Function Theorem to $±(A, 7) = 0, where $±(A, 7) = /(A) ± jF(X) 
and ^$±(z°,0) ^ we obtain a unique solution z 7 ' 1 * 1 = u = u( 7 ), 7 £ (— 70, 70), n(0) = z° 
of the equation $(A, 7) = 0, such that $(w( 7 ), 7 ) = 0,7 £ (—70,70) for some 70 > 0. 
ii) Consider the equation = -D+(A) = -D+(A) — ^s + (A)s_(A). Using 

Z^(A) = 4sin 2 ^sin 2 ^-, S± (A) = — sin ^ cos ^, cosn_(A 2n)1 ) ^ 1 (4.8) 

Z Z 7^-|- Z Z 
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we obtain 2 equations for the zeros of -D+(A) 

\ \ • V+ 7 2 V+ ,V- n • ^+ n /-/i 

OA, 7 = sin cos — cot — = 0, sin — = 0. 4.9 

v ' 2 47?_7/ + 2 2' 2 v ; 

sin — 

The equation 2 = has the zeroes A° 2n . Consider the first equation in (|4.9j) . Applying 

the Implicit Function Theorem to $(A, 7) = 0, where Jx^A^^O) ^ we obtain a unique 
solution u = m„(7),7 G (— 7 n , 7 n ), w n (0) = A° n i of the equation $(A, 7) = 0, such that 
$(tt n (7), 7) = 0,7 G (— 7 n ,7 n ) for some 7 n > 0. The proof for Dl(X) = is similar. ■ 
2. The perturbed operator. We consider the operator T y>v = —-^2 + q^,v where the 
potential q 7jl/ = aJ + -yv u Ji satisfies the Condition A, 7 G M is small and a > 0. We 
determine the asymptotics of the function p(A, ? 7)V ), At m (A, ? 7)W ), m = 1, 2. 
Lemma 4.3. Eac/j function p(A, g 7i ^), /i m (A, g 7) „), m = 1, 2, z/ = 1, |, |, zs analytic in C 2 . 
Moreover, uniformly on any compact in C 2 i/ie following asymptotics are fulfilled: 

p(A, q^) = pT(A) + o(l), /i m (A, g 7 ,,) = /^(A) + o(l), m = 1, 2, */ - 0. (4.10) 

Proof. The fundamental solutions $j, i7 , <£v i7 of the equation — y" + q^y = Xy, satisfying the 
conditions $j, i7 (0,A) = (<pj, )7 )'(0, A) = I 2 , (^ 7 )'(0, A) = <p,, j7 (0,A) = are the solutions of 
the integral equation 



<p v „(x,\)=<p (x,\)+ / <p (x - t, \)q ltV (t)<p Vjl (t, X)dt. (4.11) 

Jo 

The standard iterations in (|4.11j) yield 

= Yl ^JPwi ' <fn,u, 7 (x, A) = / y> (a5-t,A)g 7lI/ (t)^i )I/)7 (t,A)tft. (4.12) 
The last identity gives 

^0(^1, A) ( Yl (fo{x k+1 -x k ,X)q^ u (x k ))dx 1 dx 2 ...dx n , (4.13) 

1 <h<<n 

0<xi<x 2 <...<x n+ i=x i^rc^« 

where for matrices ai, a2, a n we define J] = a n a n _i...ai. Substituting the estimate 

||VA^o(^,A)|| e l Imv/X l a ' into dUHD we obtain || yfXip^^x, A) || < ( 2 »(a+l7l)) n e [im^[ a!> which 
shows that for each x ^ the formal series (|4.12|) converges uniformly on bounded subsets of 
C. Each term of this series is an entire function. Hence the sum is an entire function. Since 
v v — > p in the sense of distributions we obtain <p n;Ua (x, A) — > (p ni o >7 (x, A) as v — > uniformly 
on bounded subset of Ix C 2 , which yields (14.10}) . ■ 
We give 

Proof of Proposition 1.3. Lemma 4.3 yields p(X,q ltU ) — > p(A,g 7 , ) and p m (A, q ljU ) — > 
A*m(A, q 1 ), m = 1,2, uniformly on any compact set in C 2 as v — > 0. Then their zeros converge 
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to the corresponding zeros at v — 0, uniformly on any compact in C 2 . Due to Lemma 3.3 
we have convergence of D±(\,q JU ) and of the Lyapunov function A(A,g 7i ^). Thus using 
Lemmas 4.1-4.2 we obtain the proof of Proposition 1.3. ■ 
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